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Abstract. We siiow iiow tlie measurement of appropriately constructed particle-energy /momentum corre- 
lations allows access to the bulk viscosity of strongly interacting hadron matter in heavy ion collisions. 
This measurement can be performed by the LHC and RHIC experiments in events with high-particle 
multiplicity, following up on existing estimates of the shear viscosity based on elliptic flow. 
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1 Introduction 



A relatively broad consensus arose in the last years around 
the nature of the medium created in heavy ion collisions 
when it has become established (to the surprise of many) 
that the medium created at energy densities beyond the 
phase transition to a quark-gluon medium is best charac- 
terized as a strongly interacting fluid. This fluid features 
a very low shear viscosity over entropy density rj/s ~ 0.1, 
that is, a liquid close to perfect where neighbooring fluid 
elements strongly couple (although some alternative pos- 
sibilities such as plasma instabilities could still allow for a 
weakly coupled medium [T]). This apparent low viscosity 
has been confirmed recently at the LHC [2]. Once clear 
that the medium is behaving as a fluid, one would like to 
explore all other transport coefficients. 

Bulk viscosity is of current interest for two reasons. 
The first is that, near the phase transition, the condition 
?7 ^ C of usual fluid mechanics, that allows one to neglect 
the bulk viscosity is probably not valid in heavy ion colli- 
sions. Indeed it has been pointed out that ry/s may have 
a minimum at the phase transition [3] although it is not 
clear that the minimum must lay at precisely this transi- 
tion in an arbitrary system, as shown by the counterex- 
ample of the linear sigma model in large N [3]. Moreover, 
near the phase transition the bulk viscosity might have a 
maximum [5115], although again in certain models such as 
the one-dimensional Gross-Neveu model the bulk viscos- 
ity is monotonously decreasing [J. 

Notwithstanding the uncertainties, the scenario of a small 
shear and a relatively large bulk viscosity remains possible 
near the quark-gluon liquid phase transition, making the 
volume viscosity more important than previously thought 
from phase-shift based hadron evaluations [5]. 



The second reason is the obvious realization that the 
bulk viscosity controls the relaxation of the trace of the 
stress-energy tensor when it separates from equilib- 
rium, and is therefore related to dilatation invariance. 
Since this is a symmetry of the classical QCD Lagrangian 
(in the chiral limit), bulk viscosity will become a probe 
of the quantum trace anomaly that erases the symmetry 
from physical amplitudes and spectral properties (such as 
the proton mass). 

The most popular way of estimating the shear viscos- 
ity, by means of elliptic flow, is by its very definition how- 
ever not useful for the measurement of the bulk viscosity 
C (we will make no distinction between bulk, volume and 
second viscosity). 

Alternative ways of accessing the second viscosity have 
been proposed [9] but, in our view, the most promising 
avenue for a timely estimate is provided by the method 
of particle momentum correlations/fluctuations, that we 
study in this article, and that has been put to use to test 
models of the collision [TU] . The analogous construction for 
the measurement of shear viscosity has been developed by 
Gavin and Abdel-Aziz [TTlfT^ and there are already pre- 
liminary measurements from STAR that indicate broad 
agreement between this method and the elliptic flow esti- 
mates [131. 



2 Fluctuations of the Stress- Energy tensor 

An active degree of freedom needs to maintain an en- 
ergy of order fc^T at thermal equilibrium characterized 
by temperature T, but this energy dissipates according to 
transport equations in the medium. Hence spontaneous 
statistical fluctuations have to occur in the system at the 
right rate to reexcitc the dissipated energy to the thermal 
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level. This is the intuitive explanation of the fluctuation- 
dissipation theorem |14| , as a necessity of energy equipar- 
tition. Thus, one expects the bulk viscosity to be propor- 
tional to the fluctuations of the stress-energy tensor, and 
particularly of its trace. 

If one divides the stress-energy tensor into an ideal and 
an out-of-equilibrium part 0^,j = T^^ -'tt^u, the latter can 
be expressed for a fluid element at rest but with velocity 
gradients not necessarily vanishing, as [15] 



Tik = (^^'^,fe + Vk.i - -^SikVl^lj - (viiSik + Uk (1) 

Too = Toi = 0. 

This formula features both bulk ^ and shear viscosity 77, 
and in addition to the hydrodynamical contribution pro- 
portional to velocity gradients, a fluctuating part of purely 
statistical origin Uk- 

To experimentally access the bulk viscosity, we want 
to exploit the connection between the two-point correla- 
tion function of this fluctuating contribution to the stress- 
energy tensor tik, and the viscosities, demonstrated also 
by Landau and Lifschitz |15j . as 

{Ukiruh)tirn{r2,t2)) = (5(ri - r2)J(ti - ta) x (2) 



2T 



V {SaSkm + SimSki) + ( C - ) ^ikSw. 



(here T is the temperature). 

To separate the bulk viscosity one needs to take the 
trace over the spatial indices so that 

{tu{ri,tMr2,t2)) ^18T C S{ri~r2)d{ti~t2) . (3) 

Thus, if an experiment in heavy-ion collisions could access 
the correlations of the traced stress-energy tensor fluctua- 
tions, this would allow access to the desired bulk viscosity 
in a transparent manner. 

The route taken by Gavin and Abdel-Aziz [T^ to ac- 
cess the shear viscosity is not too different, but somewhat 
simpler since only one element of the stress-energy tensor 
needs to be measured, for example ror in cylindrical co- 
ordinates. This allows the extraction to be possible with 
only pt correlations. Instead, the bulk viscosity requires 
reconstruction of the trace of the tensor, and even in a 
perfectly central heavy ion collision where there is no vor- 
ticity and T^^ can be neglected, one still needs to measure 
at least three correlations as will be shown shortly. This 
makes the measurement more challenging, but closed for- 
mulae can be found without more conceptual difficulty. 

Let us first consider the case of a fluid element near 
rest so we can ignore relativistic corrections, and adopt 
cylindrical coordinates x' = (r, (/), z) with the z direction 
along the collider beam axis. The Cartesian tensor trans- 
forms under this coordinate change x — > x' as 



dx * dx ^ 
■J- ''J —- 

dx*^ dx'' 
The trace of the tensor is then 



^yy 



(4) 



(5) 



in cylindrical coordinates. Now, if we assume that there 
is cylindrical symmetry for central collisions, one expects 
that dv/d<f> = for the velocity (or any quantity whose 
azimuthal derivative is taken), meaning that t'^'* ~ 0. 
Thus we are left with 

{Tu{ri,ti)Tii{r2,t2)) = (6) 

{{Trr + Tzz){ri,ti){Trr + ) (r2 , ^2)) 

and as advanced this requires the measurement of three 
correlation components, due to the crossed terms, which 
will lead to the necessity of full energy reconstruction. 



3 Particle momentum fluctuations 

In this section we will show how the correlation of the 
fluctuation of the stress-energy tensor can be related to di- 
rectly measurable particle momenta. Our discussion here 
closely parallels that of Gavin and Abdel-Aziz [Tnil2| . 
with the modifications appropriate for the bulk instead 
of the shear viscosity. 



3.1 Average Equilibrium Hypothesis 

In heavy-ion collisions the experiment (e.g. a head-on lead- 
lead collision at the LHC) is repeated numerous times, 
so that a large database of equally prepared systems is 
formed. The ( ) average symbol is then understood as an 
average over all the recorded central collisions (those with 
large particle multiplicity are in common practice assumed 
to be head-on and therefore, more central). 

The fundamental hypothesis underlying the analysis is 
that a state of hydrodynamic equilibrium is reached after 
the collision. This is supported by a large body of data 
from the RHIC experiments and is widely assumed to be 
a good approximation to reality. 

The analysis we present requires a slight extension of this 
assumption. We require, as in , that the average taken 
over all collisions coincides with the equilibrium state, that 
is, the dissipative part of the stress-energy tensor aver- 
ages to zero (Tij) — 0. This excludes systematic devia- 
tions from equilibrium that may affect all collisions, and 
in effect attributes deviations from equilibrium to event 
by event fluctuations. We refer to this as the "Average 
Equilibrium Hypothesis" . With this caveat in mind, we 
will obtain next in subsection [32] the relation between the 
fluctuations of the stress-energy tensor and those of par- 
ticle momenta. 

As a last comment on this hypothesis, let us note that 
it is in the spirit of the Gibbs ensemble and the ergodic 
hypothesis, where the study of many copies of the same 
system at fixed time is equivalent to the study over very 
large times of the fluctuations of one given system. 

This large time also appears in Kubo's formula for the 
bulk viscosity 



C^ihml 

9 uj->-a u) 



dt / dxe*"*([T,,(t,x),r,,(0)]) (7) 
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that makes natural to employ Eq. ([3]) as a starting point 
to access the viscosity. 

3.2 Correlation of stress-energy tensor (tt) and 
particle momenta. 

Let us concentrate on the correlation ( ) in Eq. (|6]), 
the extension to include the zz — zz correlation and the 
rr — zz cross-correlation being trivial. 

The number of particles per unit of phase space is, in 
terms of the distribution function, 



dn = fdx 



(2^ 



(/cq + 



d^p 
(2^ 



(8) 



where we split the off-equilibrium part of the distribution 
function Sf. 

It is easiest to anticipate the answer and consider the 
following particle-momentum correlation 



(E 

all ij 



E,. 



(9) 



In this formula, the sum over i,j extends over all pairs 
of particles in a given collision event (this includes the 
square of the function for each particle). After this sum is 
performed for each individual event, the average is taken 
over all the events of the data sample. Employing Eq. ^ 
we obtain 



Crr = / dnidn2 



.{Prlf {Pr2? 



^1 



E. 



d^Pl d^P2 ,3^ ^(Prl)'^ {Pr2f 



(27r)3 (27r)3 
d^pi d^p2 
(27r)3 (27r)3 

{Prlf {Pr2f 

El E2 



d'^xid'^X2{ 



d^xid^X2 



E. 



E2 



■/(I, 2)) 



(/cql+<5/l)(/eq2 + d72)> (10) 



where the usual Boltzmann approximation of kinetic the- 
ory (valid for moderate densities) entails a factorization 
of the two-particle distribution function /(I, 2) = /i x /2. 

Additionally invoking the Average Equilibrium Hy- 
pothesis of subsection 13.11 we can ignore the terms linear 
in 5f since {5f) = (no systematic out-of-equilibrium ef- 
fects). The term independent of 5f and the quadratic term 
contribute respectively the first and second lines of 



-I- / d^Xid^X2{Trr{y^l)Trr{'X^2)) 



(11) 



with N the total particle multiplicity in an individual 
event (preferably a large number, up to several thousand 
particles at the LHC, to select central collisions). 

This equation relates the fluctuations of the stress- 
energy tensor to the particle momenta as measured in 



a detector, but it includes both the stochastic force tij 
whose correlator reveals the viscosity in Eq. ([3|), and the 
hydrodynamic part t^^"^ ^ the first terms in Eq. ([1]). One 
can substitute directly for tij in the terms linear in t 

(r,.(x)) = (t„(x)) = ( / -^M^/) = (12) 



(27r)3 E 

under the Average Equilibrium Hypothesis Q. 



3.3 Separation of the stochastic stress-energy 
contribution tu 

This leaves the quadratic terms (r/'/'^(xi,ti)T/'«''(x2,t2)) 
and (irr(xi, ii)trr(x2, ^2))- Gavin and Abdel-Aziz have been 
interested in the correlator of the hydrodynamic part, that 
satisfies a diffusion equation whose exponential solution 
has a characteristic diffusion time from which they can 
read off the shear viscosity, if one studies the long-time 
evolution of the system. In this article we also propose 
to study the correlator of the stochastic force tij whose 
intensity also gives access to the viscosity. The key to sep- 
arate them is to perform a time-integration. Since Eq. ^ 
contains a 5[ti — t2) (this just means that the stochastic 
force is statistically independent at each instant) whereas 
the hydrodynamic part features a soft time dependence 
g-At |-]^g time correlation function, an integration over 
very small time AT picks up the stochastic part. 

Indeed, for a general fluctuating magnitude satisfying 



dx 

with a continuous time-correlation function 

{x{t)x{t')) = (x2)e-^l*-*'l , 
the double integration 



(13) 



(14) 



dti I dt2{x{t)x{t')) 
AT J AT 



(15) 

is a second order infinitesimal. 

However integrating over the stochastic part on the right 
hand side of Eq. ^ (equivalent to the correlator {y{t)y{t')) 
yields 



/ dti / dt25{ti - t2) = AT 

J AT J AT 



(16) 



I AT J AT 

which is of first order. 

Thus, to separate the statistically fluctuating and hydro- 
dynamic parts we employ an integration over a short time 
interval. Under this integration one can exchange (tt) — !• 
(tt) in Eq. (HH. 



^ As an aside, and since (rrr) = 0, one has J irrrTrr) = 
/ [{TrrTrr) — {Trr){Trr))- Therefore One can equally talk of 
ATrr if any conceptual advantage could be obtained. 
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One should note that the obtained fluctuations arc al- 
ways a lower bound since part of the particles will not be 
reconstructed due to detector inefhciencies, so the actual 
bulk viscosity will be larger than the value so obtained. 
Therefore, Eq. ([T5)) and following should be understood 
as > inequalities when estimated with experimental data. 

Equation psp can be given an alternative form in terms 
of each particle's energy and mass by noting that 
(pI + p1) = - m'^ as 



(E 

all ij 



Ei Ej E 



18T C AV/AT : 



4 Relativistic boost of each fluid element 

The analysis of correlations presented hinges on Eq. (Il} 
that is valid in the fluid's rest frame. However, the fluid 
elements in the nuclear explosion following the collision 
are boosted in the laboratory frame. If the four-velocity 
of the fluid element is known to be If^, Equation (fT9)) can 
be taken to the laboratory frame by introducing the time- 
dilation factor 7 and noting that Ei = pi-U. The result 
is 



18T C j^AVi^,,/A% 



lab 



(E 



{{p-U)^-m%{{p-U)'-m%^ ^^^^{p.Uf-m\ 



iP-UUp-U), p-U 

To turn this formula into an experimentally useful ex- 
pression we need to make a few more remarks. Let us as- 
sume that one has identified a set of kinematic cuts that 
select a swarm composed of those particles coming from 
the fluid element ZiMab (this will be addressed shortly) 
during the time interval ATuh- 

The fluid element's rest frame will coincide with the center 
of mass frame. Therefore its velocity can be obtained from 
the particle swarm's energy-momentum in the laboratory 
frame as 

/3 = |^. (21) 



Then 7 = (^1 - /32)-i and 

C/ = 7(l,/3). (22) 



4 

3.4 Experimental observable in the fluid's rest-frame 

Finally we have the sought after relation between the cor- 
relations of the stress-energy tensor fluctuations and an 
observable in terms of particle momenta. 

/ / dtidt2 d^Xi(fx2{trr{xuti)trr{x2,t2)} {17) 
J AT J AT J 

\ all ij ^ ) 

where the average on the last term is a double average 
over both the event and the sample of events under certain 
kinematic cuts, multiplied by the average multiplicity, and 
the particles are supposed to have been emitted during the 
small interval AT (how this can be seen experimentally is 
left for the next Sec. E]). 

Equation (|17p is an experimental observable that can 
be constructed by measuring momenta and energies alone. 
The integral cannot be extended to the entire collision vol- 
ume, since different fluid elements have wildly different ve- 
locities, and the analysis we have performed assumes the 
local rest frame of the fluid. Since the measurement is per- 
formed in the laboratory frame, we will lift this restriction 
in Sec. El 

For the time being let us take a small clement of fluid 
Sav '^^^ characterized by a small rapidity and transverse 
velocity so that the non-relativistic analysis is a reasonable 
starting point. Then integrating Eq. ([3]) over this volume 
and the time duration of the particle emission AT, we 
have 

18T C AV/AT = (18) 

(E%^%^)-w'(S)' 

all * ^ 



^ — ^ Ei Ej E 

all ij •' 

+2(E^%^)-2W'(f)(§)- 

all ij ^ ^ 

Under the assumption of purely radial transverse flow (no 
vorticity) that we have invoked one can identify pr = p± , 
the perpendicular particle momentum. 

The right hand side of equation ([T5|) is an observable 
constructed with the momenta of all particles in the event 
database, the left hand side features the volume viscos- 
ity, that can be thus estimated if the emission time and 
volume (a hydrodynamic problem) and the temperature 
(that can be obtained by other observations such as pho- 
ton or particle spectra) are known. A theoretical observa- 
tion of interest is that the event average corresponds to 
a microcanonical ensemble average, as the energy in the 
heavy ion collision is fixed while the temperature can fiuc- 
tuate. However as in the next section we will keep a small 
volume element, we will consider its temperature as well 
defined, with the rest of the fluid acting as the heat bath. 



Once U corresponding to the fluid clement has been so 
constructed, one can compute all the products pt ■ U in 
Eq. (1201) as 

p,-U^j{E,-p,-l3) . (23) 
The four- velocity [/^ satisflcs U-U ~ 1, can be parametrized 

as 

U^ = ^y^T^Tii^coshT;, li^ cos 0, Mj^ sin0, y^r~M^sinh77^ . 

(24) 



Dobado, Llanes-Estrada and Torres-Rincon: Bulk viscosity and momentum correlations 



5 



Now let us address the fluid element's (average) space 
and time sizes AVinb, ^Tiab in terms of observable quan- 
tities. This requires understanding of the hydrodynam- 
ics of the expanding fireball, and here we will contempt 
ourselves with the simplest of models, a spherical expan- 
sion characterized by a freeze-out surface at time tj (this 
is a valid approximation if the formation radius is much 
smaller, tq ^ t/, else the polar caps of the sphere are dis- 
torted, and if the elliptic flow is moderately small). The to- 
tal swarm's longitudinal momentum will be Pz in the direc- 
tion of the heavy ion beam and is usually traded for the ra- 
pidity and the fluid element's total energy = E tanhvy, 
or 



1, f E + P, 



For our argument in this section we will consider pure ra- 
dial flow, so that the swarm's perpendicular momentum 
Pj^ in the transverse plane is parallel to p, the radial vec- 
tor in cylindrical coordinates, so the radial direction is 
automatically determined by the measurement of Pi for 
the swarm. 

We would like to express AViah and ZiTiab in terms 
of the momentum spread of the chosen particle swarm, 
centered around energy E, transverse momentum P±, az- 
imuth (j) and rapidity rj. 

Let us start by ZiTiab- We note that in the time of 
freeze-out the particle travelled a distance p — TfP± from 
the origin {f3± — P±/E). A particle arriving at the freeze- 
out distance a time ZiTiab later will have lagged by pAl3±. 
Therefore we have 



and differentiating E = mj = mj \J\ — 

_ Tf AErri^ 
ATui,^^ — — . 



(25) 



(26) 



Turning now to the spatial cylindrical coordinates, 

AViab = AzpApA(j) . (27) 

The longitudinal velocity gives Az = TfAjSz- Likewise, 
Ap = TfAl3j_. Altogether, employing again the defini- 
tion of /3 in terms of the total energy and momentum in 
Eq. (EH), 



AVi 



lab 



E 



(28) 

Finally, eliminating AP^ in terms of AP±andAE ^ we find 



AV. 



lab 



A% 



lab 



AE m? 



1 



tanh f] 



El 

E 



AE 



(29) 



tanh 1] E^ 



AP, 



^ - ^AE 
E E^ 



Note that differentiating the invariant mass of the swarm 

the three cuts AE, AP± and Arj are not independent, 
satisfying the constraint 

EAE = tanh rjE^Ai] + cosh^ rjP^AP^ . (30) 



5 Kinematic cuts 

In this section we discuss the options for the kinematic 
cuts, particularly APi_, A(f) that are workable for an ex- 
perimental collaboration, considering especially the AL- 
ICE experiment at the LHC. 

In devising them, we have to compromise between sev- 
eral constraints. 

— First, since our method calls for the separation of an 
interval AT smaller than the lifetime of the entire col- 
lision, to isolate the fluctuations, we need to consider 
a fluid element that is actually in motion and provides 
us with a clock as in section |4l Therefore we will need 
to impose a Pj_ cut that excludes Pj_ = 0. 

— Second, not all particles in a swarm move parallel enough 
to the average velocity U and may end up in a different 
element of phase space. To quantify the theory error 
introduced by this effect we have written a small Monte 
Carlo programme described shortly. 

— Third, the phase space element chosen for the mea- 
surement needs to contain enough particles across the 
collision data base to make a measurement possible. 

— Fourth and last, we have to consider that ALICE'S 
rapidity acceptance is limited to the interval (—1, 1) 
(the barrel spans about 46 degrees in polar angle to 
each side of the collision point). 

The crux of the matter is in the second point. The 
pion emission due to the freeze out of a fluid element at 
rest can approximately be described by a Bose-Einstein 
distribution in momentum p. 



dN 



= c- 



p^dp 



1 



(31) 



characterized by a temperature T and chemical potential 
p. This emission is isotropic in the rest-frame of the fluid, 
but if the fluid element is boosted, the boost velocity has 
to be compounded with the particle velocity (according 
to the special-relativistic velocity transformation rule). If 
the boost velocity is large enough, it dominates the com- 
position. Most particles are emitted aligned with /3. 
However if the boost velocity is of order of the Bose- 
Einstein velocity allowed by this distribution, the emis- 
sion becomes less beamed and each element of phase space 
is populated by particles emitted from different fluid ele- 
ments. 

In view of our fourth point above, since the longitu- 
dinal boost accepted by the ALICE detector has at most 
I77I ~ 0.9; we will consider the central part of the collision, 
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that is, take the entire longitudinal acceptance as one bin 
with 77 = 0, Arj ~ 1.8. Neglect of longitudinal momen- 
tum allows to write Eq. pip in terms of the transverse 
momentum alone as 



dN 



c- 



PldPi 



Xy/Pl+'m^-t^)/{kBT) _ ^ 



(32) 



To assess the kinematic cuts we proceed by writing 
a Monte Carlo programme. Employing Von Neumann's 
rejection method we generate a sample of several thou- 
sands of pions (corresponding to a few simulated collision 
events) distributed at random in (/) and according to the 
ALICE experimental P± distribution (TB] in 900 GeV pp 
collisions, that is well fit by and ad-hoc formula 



1 



1 d^Nr, 



11.47 e-4 ioPi P_i < 1.7 GeV 
0.25 Pr^-^^ P\ > 1.7 GeV 



(33) 

This we call defining sample and is only used to construct 
average boost velociticso 

To explore pairs of {AP± , A(j)) cuts we select the pions 
from the defining sample whose momenta fall within the 
so chosen fluid cell. We sum their momenta and energy to 
construct the cell's velocity according to Eq. (j2ip . 

Once the fluid cell has been defined and the average 
velocity is known, we turn to Eq. (j32p and generate a 
second sample of thermally distributed pions in the rest 
frame, also by Von Neumann's rejection method, the ther- 
mal sample. 

This sample represents isotropic emission in the fluid's 
rest frame and we impose no restriction on P±^ or ex- 
cept thermal distribution. 

Finally, we apply the Lorentz boost with the velocity from 
Eq. (PT|) corresponding to the fluid cell to each of the pions 
in the thermal sample^ and examine what fraction of them 
falls outside of the initial kinematic cuts that defined the 
fluid cell. 

We find that a non-negligible but controllable percent- 
age of the sample pions end up into a different fluid cell. 
The results are listed in Table [1] as percentages of particles 
appearing with momenta that would correspond to a fiuid 
cell other than used to generate them. 

For completeness we also address ALICE'S Pb+Pb data 
at y/s = 2.76 TcV. We fit the P± distrbution in analogy 
with Eq. by 



1 



1 d^Nch 



Ne^ut 27rF_L d7jdP± 



8.90 • 10^ e-2-93P_L < 2.0 GeV 
2.00 • 10^ p-6-29 > 2.0 GeV 

(34) 

to obtain the corresponding defining sample and repeat 
the analysis (obtain each cell's velocity, generate a thermal 
sample, boost the pions thereof and examine their final 
momenta) . The corresponding result is given in table [21 
Examination of table [T] teaches several general lessons. 



^ Incidently, the same data taken at low P± can be used 
to fit the rest-frame Bose-Einstein thermal distribution param- 
eters (temperature and pion chemical potential) in Eq. p2|) . 



— If the boost velocity is generally larger (the average 
momentum is at higher P_l), pions do not spread out 
too much and losses from the cell are lowered. 

— If the azimuthal-angle cut A(j> is larger, losses from the 
cell are in general smaller because, after boosting the 
thermal sample, most pions remain inside this larger 
cone. 

— If on the other hand the azimuthal-angle cut is very 
small, low momentum particles find it easy to leave the 
tiny resulting angular cone. One can reduce the mixing 
between fiuid cells by proceeding to larger P± so the 
boost focuses the swarm in the correct direction. 

— In the extreme case, if the momentum cut is centered 
at huge momenta, the cell's /3 is very close to 1. Al- 
most independently of the initial thermal configuration 
most of the particles follow the boost and fall within 
the defining momentum cut. By increasing the angular 
acceptance this proportion is further improved. How- 
ever the statistics with real data fall exponentially with 
P±, so a balance has to be found between larger mo- 
mentum and sufficient data. (At too large momentum 
one should not trust thermalization either). 

A reasonable choice would be for instance to take a 
small angular cut of 60° and identify all pions with P± > 
0.3 GeV. The number of particles that mix with other fluid 
cells is then around a third. This mixing should be consid- 
ered a systematic theory uncertainty in the measurement 
of the bulk viscosity. 



6 Conclusion 

We believe we have identified a way that allows to access 
bulk viscosity from energy-momentum correlations. The 
modus operandi that we suggest would be to define three 
appropriate kinematic cuts A(j), AP± and AE defining a 
swarm of particles centered around (p, P± and P to a set of 
recorded central collision events (many such swarms can 
be defined and the results compared). 

Substituting then equation ([^^ in ([^0)1 and evaluating 
the right-hand side correlator over the data, one obtains 
an estimate for the bulk viscosity 



E'^AEm^ , /(p-U) 



18Tj^TjA(j)Pl 
1 



p-U 



{id-V, - P^) - ^AP^] {EAP^ P^AE) 

(35) 

that depends on the temperature T and the freeze-out 
time Tf . These can be obtained from other measurements 
and then grant access to the bulk viscosity. By varying the 
size of the cuts AE, Ac/), AP± one can explore the attend- 
ing systematic uncertainties, and by varying the central 
value of these variables around which the particle swarm 
is chosen, one can study the variations of the bulk viscosity 
over the collision volume. 
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Table 1. We show average velocity /? of the given swarms of particles within the azimuthal angular Z\(j!> and the transverse 
momentum AP± kinematic cuts, with the particles distributed according to Eq. H33[l corresponding to proton-proton collisions. 
We also show, for each given binning with velocity /3, the percentage of thermally emitted particles following Eq. (|32|) that are 
lost from the bin after compounding /3 with the particle's thermal velocity. Typical results show that a fourth to a third of 
particles with well-chosen cuts populate other fluid elements introducing an irreducible theory error. 





All Pj_ 


P± > 0.3 GeV 


P± > 0.5 GeV 


P± e (0.3, 2) GeV 


P± e (0.3, 3) GeV 




/3 % 


P % 


/3 % 


/3 % 


/3 % 


Acf) = ±20° 


0.93 41.6 


0.96 36.5 


0.97 36.7 


0.96 62.0 


0.96 48.8 


Zi0 = ±30° 


0.91 31.9 


0.93 33.3 


0.94 36.9 


0.93 51.0 


0.93 39.9 


A(t> = ±45° 


0.86 24.7 


0.88 33.7 


0.89 42.4 


0.88 42.2 


0.88 35.6 


A(t> = ±60° 


0.79 18.6 


0.81 35.7 


0.82 49.2 


0.81 39.6 


0.81 36.3 



Table 2. Same as in table [T] but for lead- lead collisions, with the pion distribution following Eq. (|34[) . 





All Pj_ 


Px > 0.3 GeV 


Pj_ > 0.5 GeV 


Px G (0.3, 2) GeV 


Px e (0.3, 3) GeV 




P % 


P % 


P % 


P % 


P % 


A(f> = ±20° 


0.95 38.1 


0.96 36.1 


0.97 36.6 


0.96 69.9 


0.96 57.2 


A(t> = ±30° 


0.92 30.0 


0.94 33.1 


0.94 36.4 


0.94 59.7 


0.94 47.3 


All) = ±45° 


0.87 23.9 


0.88 33.2 


0.89 40.9 


0.88 49.9 


0.88 39.7 


A(j) = ±60° 


0.80 18.2 


0.81 34.3 


0.81 45.8 


0.81 44.1 


0.81 36.8 



Alternatively, for midrapidity (r/ ~ 0) one can use the 2. 
approximate formula: 

3. 



c 



1 



1 



18T-f^T^ A(f)AT] Pi - Pi 



A 



p-U 



,(36) 



7 — f— , - X " - ± \ ^ - / 4 

where for the ALICE detector the longitudinal acceptance 
is At] ~ 1.8. Importantly, this factor is independent of 
APa_. 5 

One may wonder whether the bulk viscosity measured 
from Eq. psp carries information about the entire colli- g 
sions process or only the late stage near kinetic freeze- 
out. Since we have chosen ZiTiab to pick-up fluctuations j 
in Eq. , the measurement refers to precisely the inter- ^ 
val Z\7iab before freeze-out. For exmaple if AE ^ 2 GcV, 
E - 500 McV, m - 140 MeV and r/ - 10 fm/c, then 9 
ATlnh refers to the last 3 fm/c before the freeze-out. 

10, 



11. 



12. 
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